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INTEGRATED FUNCTIONALS OF NORMAL AND 
FRACTIONAL PROCESSES^ 

By Boris Buchmann and Ngai Hang Chan 
Monash University and Chinese University of Hong Kong 

Consider Z/(m) = /J" f{Ns) ds,t>0,ue [0, 1], where N = {Nt)teR 
is a normal process and / is a measurable real- valued function satisfy- 
ing Ef{No)^ < oo and Ef{No) — 0. If the dependence is sufflciently 
weak Hariz [J. Multivariate Anal. 80 (2002) 191-216] showed that 
Z{ jt}!"^ converges in distribution to a multiple of standard Brown- 
ian motion as t ^ oa. If the dependence is sufficiently strong, then 
ZtjiEZtiXf'Y^'^ converges in distribution to a higher order Hermite 
process as t — > oo by a result by Taqqu [ Wahrsch. Verw. Gebiete 
50 (1979) 53-83]. When passing from weak to strong dependence, a 
unique situation encompassed by neither results is encountered. In 
this paper, we investigate this situation in detail and show that the 
limiting process is still a Brownian motion, but a nonstandard norm- 
ing is required. We apply our result to some functionals of fractional 
Brownian motion which arise in time series. For all Hurst indices 
H G (0,1), we give their limiting distributions. In this context, we 
show that the known results are only applicable to H < 3/4 and 
H > 3/4, respectively, whereas our result covers H = 3/4. 

1. Introduction. With the increase popularity of using self-similar pro- 
cesses to model econometric time series (see Bailie [2]) and network traffic 
(see Willinger, Taqqu, Sherman and Wilson [27]), understanding the limit- 
ing behavior of integrated functionals of normal processes becomes a chal- 
lenging and important task. Functional limit theorems of this nature have 
found applications in various disciplines (see, e.g., Bhattachary, Gupta and 
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Waymire [4] for an interesting application in hydrology). One of the early 
results of this nature is given in the seminal paper by Taqqu [25], where it 
was shown that the limit of such processes are in the domain of attraction 
of higher order Hermite processes when the underlying normal process ex- 
hibits strong dependence. This problem was subsequently studied in Hariz 
[15] when the underlying normal process is weakly dependent. The problem 
becomes more intriguing when transitions from weak to strong dependence 
take place. Among others, we establish the limiting distribution of integrated 
functionals of fractional Brownian motions for all Hurst exponents in (0,1) 
in this paper. 

Let {Q, P) be a probability space carrying a standard Brownian motion 
W = (Wj)tgK and let a be a Borel function with / a^(s) ds = 1. We consider 
a normal process N defined by 

(1.1) Nt= [ a{t-s)dWs, t£M, 

and also its covariance function r : R ^ M, where we set 

r{t) = E{NtNo), t>0. 

For /c G No, let 7ik denote the kth Hermite polynomial with leading coeffi- 
cient one defined by 

(1.2) = (-If exp(f)-^exp(4). 

In the sequel, let / be a measurable real-valued function such that / satisfies 
EfiNo) = and Ef'^^No) < oo. Then / can be expanded into an orthogonal 
series as 

oo 

(1-3) /(x) = ^|Wfc(x), 

k=l 

where the coefficients are given by 

(1.4) Ck = Enk{No)f{No), ken. 

The series in (1.3) converges in Lp'{v) where v denotes the standard normal 
law on the Borel sets of M. In this case, the constant 

g = g(/) = inf{fc:cfc/0} 

is called the Hermite rank of / by Taqqu in his pioneering paper [24]. 

The Hermite rank plays a crucial role in the asymptotic behavior of 
sample means of nonlinear functionals for both Gaussian processes and 
fields indexed by continuous and discrete parameters (cf. Breuer and Ma- 
jor [7], Chambers and Slud [11], Dobrushin and Major [13], Hariz [15] and 
Taqqu [24] and [25] and references therein). 
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In particular, Hariz [15] showed that, if 

POO 

(1.5) / |r(s)|«ds < oo, 

Jo 

then a central limit theorem (CLT) holds, that is, the finite dimensional 
distributions of zl{-)/t^l^ = Jq /(iV„) d-u/i^a 

converges to aW, as t — > oo, 

where 




The integral in (1.5) is finite if the covariance function of N tends to zero 
sufficiently fast, that is, the dependence of the underlying process N is suffi- 
ciently weak. When the underlying process is strongly dependent, a noncen- 
tral limit theorem (NCLT) holds: Under some technical conditions on the 
kernel a, Z-^ is in the domain of attraction of a Hermite process of order 
q. This was shown by Taqqu [25]. However, there are situations where the 
quantity in (1.5) is infinite, but the dependence is still too weak to satisfy a 
NCLT. This situation typically occurs when one wants to explore the domain 
of attraction of integrated functionals of fractional Ornstein-Uhlenbeck pro- 
cesses for all Hurst indices H S (0, 1). As a standard rule, the CLT covers 
H £{0,1- l/(2g)) whereas the NCLT is applicable to H £ {I - l/{2q), 1) 
[q> 2]. In this paper, we are interested in the boundary corresponding to 
H=l-l/i2q). 

The paper is organized as follows. In Section 2, we state a general theorem 
dealing with the boundary case. Also, we give an application of our result 
to fractional Ornstein-Uhlenbeck processes. An important application of 
results of this nature arises from the study of strongly dependent times 
series (cf. Buchmann and Chan [8]). In Section 3, we review and extend 
the asymptotic theory of the so-called unit-root problem as an application 
of Section 2. In this context, the CLT and the NCLT are only applicable 
to H £ (0,3/4) and H £ (3/4,1), respectively, whereas our theorem covers 
H = 3/4. All proofs are found in Section 4. 

2. Main result. Throughout let C(]Rq') be the space of continuous func- 
tions on Wq = [0,oo) endowed with the convergence of locally uniform con- 
vergence. The proof of the following theorem is found in Section 4.1. 



Theorem 2.1. Suppose that (1.1) is satisfied. Let L and be slowly 
varying at infinity with limj^oo L{t) = oo and 

L\At) 

(2.1) limsup ' ' < oo, 

t^oo L(t) 
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where 

(2.2) L{t)= f r'i{u)du, UM = I \r{u)\'i du, t>0. 

Jo Jo 

Let f £ L'^i'i^) with E/^Nq) = and Hermite rank q = q{f) > 1. Then the 
following assertions hold: 

(i) In the sense of the finite dimensional distributions, 

(2.3) (2tL(t))- V2 /(iv„) du 4 -^W, t ^ oo. 



(ii) In addition, if there exists R> 1 such that J2'h=q ItIi-^^ finite then 
(2.3) holds in C(M(j"). 

Remark 1. (i) It follows from Karamata's theorem, Fubini's theorem 
and the diagram formula (A. 2) that 



E 



f ng{Ns)ds =2q\ I I r'^{u)dudv r^2q\tL{t), t ^ oo, 
Jo J "'0 "'0 



and thus L is eventually positive at infinity. 

(ii) To set up weak convergence to a Brownian motion as in (2.3), it is well 
known that a possible norming function has to be regularly varying with in- 
dex 1/2 (cf. Lamperti [17] and the discussion in Taqqu [24]). Consequently, 
the assumption that L is itself slowly varying is rather weak. Hariz [15] gives 
a functional version of his result under similar conditions on the Hermite 
expansion of / and A^. The existence of a representation of form (1.1) is 
equivalent to the absolute continuity of the spectral measure of N with re- 
spect to Lebesgue measure (cf. Ibragimov and Linnik [16], Theorem 16.7.2). 
This covers a lot of important cases and, in particular, Ornstein-Uhenbeck 
processes driven by the fractional Brownian motion. Our proof of Theorem 1 
relies on elegant results obtained by Peccati and Tudor [20] and formulated 
in terms of functionals on a Gaussian space. Therefore, we need this space 
controlled by a diffuse measure, which is taken to be the Lebesgue measure 
in the proof. It would be interesting to investigate whether the arguments by 
Peccati and Tudor [20] extend to other measures than the Lebesgue measure. 
However, this is beyond the scope of this paper. 

(iii) For normal processes with a bounded spectral density. Chambers and 
Slud [11] have shown central limit theorems for more general functionals, in- 
cluding those depending on infinitely many coordinates. For the long-range 
dependence situation considered in our paper, the spectral density is usu- 
ally unbounded, however. The result presented in this paper only covers 
functionals depending on one coordinate. A general theory has yet to be 
developed for general functionals of long-range dependent processes. 
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To prepare our analysis in the next section, we conclude this section with 
an application to fractional Ornstein-Uhlenbeck processes. Let H G (0,1). 
A fractional Brownian motion (FBM) = {BP)t<z^ with Hurst index H 
is a centered Gaussian process with almost surely locally Holder continuous 
paths of any order strictly smaller than H and covariance function 

(2.4) EB^B^ = i(|t|2^ + -\t- 5|2^), s,t G M. 

A choice H = 1/2 relates to a standard Brownian motion W (cf. Bhat- 
tacharya and Waymire [3] for interpretations of Hurst index and Samorod- 
nitsky and Taqqu [23] for further properties of FBM) . 

In particular, FBM is self-similar, that is, for all c G M, in the sense of 
finite dimensional distributions, 

(2.5) {B^,)^\cf{B^). 

Also, let Bi^ = {Bi^i)t>Q be the fractional Ornstein-Uhlenbeck process (FOUP), 
defined by 

(2.6) B^|^ = e-^(*-") dB^ , i > 0, 7 G M, 

where the integral in (2.6) converges pathwisely in the usual Riemann- 
Stieltjes sense (cf. Buchmann and Kliippelberg [9] and Cheridito, Kawaguchi 
and Maejima [12] for further properties of FOUP). 

For 7 > 0, the process B^^ is asymptotically stationary. In the next corol- 
lary, we apply our result to the case when the Hermite rank q = 2. Up to 
some technicalities (cf. Section 4.2), it follows from the CLT by Hariz [15] 
and the NCLT by Taqqu [25] for H G (0,3/4) and H G (3/4, 1), respectively. 
Theorem 2.1 deals with the remaining case H = 3/4. We conclude this sec- 
tion with a remark where we summarize the corresponding results for the 
general case and allow for arbitrary Hermite ranks q. We refer to Taqqu [25] 
and references therein for the definition and properties of Rosenblatt pro- 
cesses with Hurst index H. Throughout F denotes Euler's Gamma function. 



Corollary 2.1. For all H G (0, 1) and 7 > 0, it holds in C{M.q) 
rt- r 



(2.7) gnit) 







H .2_r(2^>i) 



ds hH{'y)crHX^ , t—>-oo, 



where, for H G (0, 3/4], is a Brownian motion and for H G (3/4, 1), 
is a Rosenblatt process with E^Xi)"^ = 1 and Hurst index H . 
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The quantities gn, cth nn-d hfji'j) are given by the following formulas: 

t-^l\ 



(2.8) gnit) 

(2.9) huil) 



{t\ogt)-^/\ 

1-2H 



t 



'\I2-2H 



if G (0,3/4), 
if = 3/4, 
if G (3/4,1), 

(0,3/4), 
i/G [3/4,1), 



t > 1, 



(2.10) 



'(2/7r)^/2p(2/7+l 
X sin(7rff) 



t2-m 



(l+a^ 
3/4, 

U[(4ii-2)/(4F-3)]i/2, 



1/2 



i/G (0,3/4), 

^ = 3/4, 
if G (3/4,1). 



Remark 2. Let g G N, if G (0, 1) and 7 > 0. In Corollary 2.1, we only 
deal with Hermite rank g = 2, which is sufficient for Section 3. It is natural 
to ask whether Corollary 2.1 can be extended to other Hermite polynomials. 
The answer is affirmative and as the arguments are very similar to the proof 
of Corollary 2.1, we only state the main results here without proof. To this 
end, let 



(2.11) 



lH \2 _^ „,-2H , 



He (0,1) 



and observe that E{Bi^f)'^ —>■ IJ'Jj^ (cf. Buchmann and Kliippelberg [9]). 

In view of Lemma 4.1(iii) below, it follows from Taqqu [25] that, for q>l 
and > 1 — ^, as t ^ 00, in C(]Ro"), 

,,(i-H)-i f\^i^H^fBi^^,)ds 



(2g!)i/2[(2// - 2)q + 1] [{2H - 2)q + 2]" 



-1/2 



2H-1 
T{2H) 



q/2 



where B*^'^ is the Hermite process of order q with Hurst index H and unit 
variance E{Bf^)'^ = 1 (cf. [25] for the definition of higher order Hermite 
processes). Note that B^'^ = B^ for q = l and H £ (1/2, 1) whereas, for q = 
2 and G (3/4, 1), B^'^ = is the Rosenblatt process in Corollary 2.1. 

If g > 2 and H = 1 — l/{2q), then it follows from Theorem 2.1 and Lemma 4.1(iii) 
that, as t — > 00, in CfMo"), 



(tlogt) 



-1/2 



'n,{^^f.fBi^Jds^{2ql)y' 



2H-1 
T{2H) 



q/2 
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If either, both q = l and H < 1/2, or, g > 2 and H <1 — l/{2q), then we get 
from Hariz [15] and Lemma 4.1(iii) that, as t ^ oo, in C(]Rq'), 



(2.12) ^1/2 
where 



(2(7!)i/2 



[2r(2i? + l)]'?/2 V 7 



q,H 



1/2 



T{2H + l)e-* 
+ 2H 



dt. 



Here, the integral converges absolutely, provided either q = 1 and < H < 
1/2, or g > 2 and < < 1 - l/(2g). Observe that = 2'i/q for all 

g G N. A more refined analysis shows that Iq^n > for all g' > 2 and < H < 
l-l/(2g). 

However, note that Iih = for all < H < 1/2 such that the limit in 
(2.12) is trivial in this case. To find a nontrivial limit, consider a function 
V': [0,oo) ^ M with -0(0) = 0, that is, locally Holder of order /3 > 0. For ah 
r > there, thus exists Ct G K such that Itl^iu) — il^{v)\ < Ct\u — for all 
< u,v < T. In this case. 



sup 

0<D<T 



7 



< t' 



7 



r(/? + l)+ sup s^e~ 

0<s<oo 



for all j,t,T > 0, and thus for all 7 > 0, locally uniformly in v>0, 



(2.13) 



7 



t — > 00. 



Let H £ (0, 1), in view of (2.5), we get from simple substitutions and partial 
integrations that, for all t>0, 

(2.14) r B^,ds = t'-'' f'e^'^'-^Blids^t f e^'^'-^ ds, 



where the last identity holds in the sense of finite dimensional distributions. 
Clearly, we have Bq = with probability one. Further, sample paths of B^ 
are locally Holder of any order strictly smaller than H, almost surely. Hence, 
the right-hand side in (2.14) converges to B^ locally uniformly, almost 
surely, as t — > 00, by means of (2.13), and thus in C(Mo')) as t — > 00, 



Bi^^^ds^r'B" 
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3. The unit root problem. We apply Corollary 2.1 to some functionals 
of FBM that occur in the unit-root problem in times series analysis. To 
illustrate this problem, let e = (en)nGN be a sequence of independent and 
identically distributed random variables with variance one and consider the 
first order autoregressive model, 

Xn = PXn-i +en, n G N, Xo = 0. 

The parameter /? is unknown and has to be estimated from the observations 
Xi , . . . , Xn . The least squares estimator for /3 is given by the formula 

En— 1 V V 
n -n , ^ ^ f=0 Xt+lXt 

On - On[^0, ^n) - ^-1 ^2 ' 

Z^t=0 

For < 1 (stationary regime), Mann and Wald [18] showed that 

1/2 

(3.1) f„=IVxn {bn-P)^Wi, n^oo. 




For the explosive case |/?| > 1, (3.1) holds whenever e is a sequence of inde- 
pendent standard normal random variables (Anderson [1]). 

On the other hand, (3.1) fails to hold for /? = 1 even when e forms a 
sequence of independent standard normal random variables. White [26] and 
Rao [22] showed that it is a functional of Brownian motion W = {Wt)o<t<i, 
that is, 



(3.2) fn^fiO) = '^[W?-l]\[ 

Uo 



r-i n-1/2 
W^ds 



n — > oo. 



This contrast between (3.1) and (3.2) is a typical example of a critical phe- 
nomenon. The parameter value /3 = 1 comprises a singularity and there is a 
lack of smooth transition of the limiting distribution of f„ when /? is close to 
one. For finite sample analysis or tests under local alternatives, a key ques- 
tion becomes that if /? is close to one, what kind of approximation should 
be used for f^? An answer to this question is given in Chan and Wei [10], 
where the following class of models is proposed: 

(3.3) 4")=/3„xf!i + et, nGN, xj^^ = 0. 

Suppose that there exists 7 G M such that Pn = ^ — l/n. Then it is shown 
by Chan and Wei [10] that 



(3.4) f.4f(7) = ^" 7-^-^"% n-00. 



1/2 

where = B-y is the Ornstein-Uhlenbeck process driven by Brownian 
motion and the integral on the right-hand side is defined in the Ito's sense. 
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In [8] , a generalized asymptotic theory of ordinary least squares estimators 
for a large class of possibly strongly dependent noise sequences was given. 
The result covers certain stationary and ergodic sequences e = (e„)„gN with 
mean zero and finite variance cj^ = Ee^ , where the associated partial sum 
process is in the domain of attraction of for some H G (0, 1). In this 
case, there exist matrices D„ = Dn{H) G M^^^, n G N, G (0, 1), such that, 
as n ^ oo, 

('\^\ D ( ^ , (^H,^{l)\ . ,^VrH,i(7) 

^^•^^ Un - /3„ j ^ ^^^^/^ KthAi) ) " ^-^/'T Uh,2(7) 

where 1a denotes the indicator function of some set A. Further, TH,i{j) 
denotes the ith component of the vector th{j) G defined by 

(3.6) th(7)= (eH(7)',(^^+7^'«s)'rf«)eH(7)')', 



(3.7) 6^(7) 



where x' denotes the transpose of a vector x G M". 
It follows from Ito's formula that 

(3.8) f (7) = ri/2,3(7) - \Ti/2,i{n), 

and thus the limit in (3.4) can be derived from (3.5) and (3.8), provided 
0-2 = 1 (cf. [8] for details). 

It was shown in [10], that t(7) forms a continuous family indexed by 
7 G [—00,00], in particular, we have with (3.1) in view, 

(3.9) r(7)4p^i, |7|^oo. 

Next we want to extend (3.9) to all H G (0, 1). Therefore, we consider the 
vector Tni'y) and its asymptotic properties for I7I — > 00. In Remark 5 below, 
we show how the limit in (3.9) is derived from the general theory. 

There are two cases, that is, 7 ^ 00 or 7 ^ —00. Both cases exhibit quite 
different qualitative behavior. With (3.9) in mind, this seems to be puzzling 
at the first glance. However, the underlying dynamics are rather different 
as formally taking 7 to +00 or —00 in (3.3) corresponds to an infinitesimal 
return to the stationary regime or the explosive one. 

We first tackle 7 — > 00. The corresponding theorem is implied by Corol- 
lary 2.1 as it follows from the self-similarity of FBM that we may rewrite 
th i'y) into a functional of integrated squares of . In view of Corollary 2.1 
for Hermite rank q = 2, we thus expect two types of limit distributions de- 
pending onif<3/4or//>3/4 with an additional logarithm in the norming 
functions for H = 3/4 (cf. Section 4.3 for a proof). 
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Theorem 3.1. Let {Z,Y) a standard normal random vector. For H £ 
(3/4,1), let Rh =R^{1) be a Rosenblatt distributed random variable with 
ER^ = 1 , independent of Y . 



(i) For H e (0,3/4] and 7 > 1, there exist T.h G M.'^'', bnil) G 



anc 



DH{7)eR'^'"^ such that 



Dh{i){th{i) - bnil)) 



Z 

y2 I > 



7 00. 



(ii) For H G (3/4, 1) and 7 > 1 there exist T^h G M''^'^^ 6^(7) G M'* anc 
Dh{i) G M"^^"^ such that 



Dh{i){th{i) - buil)) 



Rh\ 

y2 ), 



7 ^ 00. 



Remark 3. For 7 > 1 and H G (0, 1), the scaling matrices Dni'y) have 
the following forms: 



diag(y/2-^,7^/^-^^,7^-^/^l), i^G (0,3/4), 

diag(7-^/^log7)"'/^7-'(log7)^'/^ 

y/4(iog^)-i/2,i), H = 3/4, 

diag(f"='^,7^-^^,V-^,l), i^G (3/4,1). 



For 7 > and H G (0, 1), the matrices T,h and centering vectors S.Hil) are 
given by the following formulas 



3/2 



/ l^H 1 \ 



\ 



KH 

1/2/ 



^1^(7) 



V 



7 



where is defined in (2.11) and we set 

roo C2-4H 



7r-^/2 sin(^//)r(2/7 + 1)^/2 
(3/8)1/2^-1/4 



-1/2 



\4:H-3 



1/2 



r(2ff- 1) 



(i+a- 



-1/2 



■d^ 



1/2 



HG (0,3/4), 
H = 3/A, 
He (3/4,1). 



We conclude this section with a corresponding result dealing with the 
other case, namely 7 ^ — 00, together with some remarks (cf. Section 4.4 for 
its proof). 

Theorem 3.2. Let He (0,1). Then for -1 ^ -00, 



diag(|7|-('^+^)/'el^l,|7|-'^-'e2H,|7|('^-^)/M7r'e''^')r//(7) 



-2H-1 2|7| I |(2H-l)/2 
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^ "HfmTT)^' iWTT) • ^('^ + ^'""•') 

X (|Z|-i,Z-2,S'5ign(Z),y/Z)', 

where {Y, Z) is standard normal and 

sign(x) = 1, X > 0, sign(x) = — 1, x < 0. 

Remark 4. In Theorem 3.2, observe that Y/Z is a Cauchy random 
variable. A similar result was obtained by Anderson [1] for the ordinary 
least squares estimator bn for nonstationary AR(1) models, that is, |/3| > 1, 
with independent standard Gaussian innovations. 

Remark 5. As indicated, we aim for recovering the limit in (3.9) by 
Theorems 3.1-3.2. An easy computation yields = 2 and ki/2 = 1/2, by 
virtue of the identity dS,/{l + = it/A. By Theorem 3.1, 

and thus by (3.8), 

r(7) = (nA3(7) - 2-^^'/') - i(n/2,i - 2V2^V2) S Z, 7-00. 

On the other hand. Theorem 3.2 yields 

(^l7|-^eHri/2,i(7)\4/ 2\Z\-^ \ ^ ^ _^ 
\ nA3(7) J [YsigniZ))^ ^ 

By means of (3.8), we get, for 7— > —00, 

t(7) = Ti/2,3(7) - ^|7|e-l^l(l7r'el^lri/2,i(7)) ^ ^ sign(Z) ^ Z. 
4. Proofs. 

4.1. Proof of Theorem 2.1. (i) First, we wish to establish (i) for the spe- 
cial choice / = Tiq, g G N. We make use of a CLT for vector- valued multiple 
stochastic integrals by Peccatti and Tudor [20] (cf. [19], Section 1.1, for the 
necessary background). Secondly, we extend the result by a reduction, this 
being similar as in [24]. 

To this end, fix q € N and let 

Ag^t{vi,V2) = {2tL{t))-^'^ ngiNu)du, V2>vi>0. 

J Hit 

In order to show (2.3) for the finite dimensional distributions for / = Tig, it 
suffices to show that, as t — > 00, 

(4.1) iAq^t{v2j-l,V2j))l<j<d l<j<di 
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for all d G N and V2d > V2d-i > V2d > ■ ■ ■ >V2> vi>Q. 

Fix d G N and V2d > 'f2d-i ^ '^2d > ■ • ■ >V2 > vi > 0. By Proposition 1.1.4 
in [19] and a stochastic version of the stochastic Fubini's theorem, it follows 
from (2.2) that, almost surely for 1 <j <d, 

A,i('y2j-l,W2j) 

(4.2) 

= {2tL{t))-^l^ / X{a{t-u,)dW^'^{u^,...,Ud), 

where the right-hand side is understood in terms of multiple Wiener inte- 
grals (note that the Hermite polynomials in (1.2) have a different leading 
coefficient than the ones in [19]). 

Suppose that we can show that, for all V4> V3>V2> vi>0, as t — > oo, 

(4.3) EAltivi,V2) ^ qlE{W,, - W,,)^ 

(4.4) E[Aq4vi,V2)Ag,t(.V3,V4)]^0 
and 

(4.5) EAl,{vi,V2) ^ q\E{W,, - W,,)\ 

Then (4.1) is immediate in view of (4.2) and Theorem 1 in [20]. To show 
(4.3), we make use of the diagram formula, which together with some nec- 
essary graph-theoretical concepts, are stated in the Appendix. 

Fix V2 > vi >0. As r is symmetric and L is slowly varying, we get from 
Fubini's theorem and Karamata's theorem that, as t — > oo, 

r'^{ui — U2) dui du2 = 2 / L{u) du ~ 2{v2 — vi)tL{t), 



I tvi J tvi 

(4.6) 

and thus by (A. 2), as t ^ cxd, 

rtV2 rtw 

EAlt{vi,V2)=2ql{2tL{t)y^ / r^iui - U2) duidu2 ^ ql{v2 - vi). 

J toi J tu 



This completes the proof of (4.3). 

Next pick f4 > fs > ^2 > f 1 > 0. By simple substitutions, observe that 

r'^{ui — U2) dui du2 



L{u) du, 



tV3 Jtvi 

t{v4-vi) rt{v3-V2) rt{v4-V2) rtivs-vi)' 



/ +/ "/ 

JO Jo Jo 



UO JO Jo Jo 

for t > 0, and thus by means of Karamata's theorem, as t — > 00, 

/'tV4 rtV2 

(4.7) / / r''{ui-U2)duidu2 = o{tL{t)), t^oo. 

J tV3 J tvi 
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Now (4.4) is implied by (A. 2) and (4.7) since, as t — > oo, 

E[Aq^t{vi,V2)Ag^tiv3,Vi)] 

= (2g!)(2tL(t))-M / r'i{ui-U2)duidu2^^. 

To show (4.5) let e > and recall (2.2). It follows from (2.1) and lim^^oo L{t) = 
oo that / := /q°° \r{u)\^ du is infinite. By assumption, L|.| is also slowly vary- 
ing. Consequently, by Holder's inequality and Karamata's theorem, for all 
v>e>0, 

j'j \r{u)r du < (rf)i-^'"/u;:;/^(et) + [{v - e)*]^-™/" [L|.| [vt) - L,., (et)]™/^ 

~ (rf)i-™/'?Lj^/^(t), t ^ oo, 1 < m < g, 
and thus by (2.1), 

(4.8) / |r(n)rd'u = o(ti-'"/«L(t)'"/''), t ^ oo, 1 < m < g, u > 0. 

JO 

On the other hand, it follows from the inequality between geometric and 
arithmetic mean and simple substitutions that for all m G N, u > and 

Q<Wi,...,Wm<Vt, 

rvt pvt 

(4.9) / W\r{u-Wi)\du<2 \r{u)\"'' du, t>0. 
Jo f-J^ Jo 

We return to the proof of (4.5). Fix V2 > vi > 0. Using the notation of the 
Appendix, let P(4,g) be the set of all diagrams with 4 levels and rows of 
length q. Let be the subset of all D £ 'D{4,q) such that the following 
implication holds for all e, / G V(D): 

m(e)=m(/) ^ M{e) = M{f). 

If D G "^(4, q) \ then, for 1 < i < 3 and i < j < 4, there exist integers 
< 

— Qi additionally satisfying mi i + ni\ 2 ~l~ ^711,3 ~l~ ^i,4 — Q and 1 < 
^712,3 + "i2,4 < <?, such that 
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Thus, by (4.8)-(4.9) and (2.1), as t ^ oo, we have 

rtV2 rtV2 rtV2 rtV2 

/ / / / n r{u^(^^)-UM{e))duidU2dU3dU4: 

JtVl JtVl JtVl JtVl 

rtV2 rtV2 

<8u2tL\.\{v2t) \r{u)r^-^dui |r(n)r2.3+m2,4^^^ 



= o(i^L2(t)), t^oo. 
Consequently, we get from (A. 2) and Fubini's theorem that, as t — > oo, 

rtv2 /•tv2 \ 2 



{vi,V2)=5{ql)H / r'^{ui-U2)duidu2) +o{t^L\t)). 



(4.10) hm / |r('u)|'^+^dn/L|.|(t) = hm |r(t)| =0. 



This completes the proof of (4.5) in view of (4.6). 

To summarize, the assertions (4.3)-(4.5) are all in place, and thus (4.1) 
is implied by (4.2) and Theorem 1 in [20]. Further, this completes the proof 
of (i) for the choice / = Jiq. To extend (2.3) to all / G recall that 

limt_»oo -^(i) =oo, and thus limt_^oo -^|.| (i) = oo by (2.1). Further, it follows 
from (2.2) that a spectral density exists (cf. Ibragimov and Linnik [16], Theo- 
rem 16.7.2). In particular, r must be continuous and we have limt^oo '"(i) = 
by means of the Riemann-Lebesgue lemma. If \r{u)\^^^ du is infinite, 
then we get from I'Hospital's rule that, for all f > 0, 

t— too Jq I - ' ' ' I I » ' t— »oo ' 

Otherwise, if /q°° |r(ii)|'?+^ du is finite, (4.10) is trivial. 

On the other hand, it follows from stationarity and Cauchy-Schwarz in- 
equality that |r(i)| < t for all t e M, and thus for all u, t > and m > g + 1, 

vt rvt rvt 

/ \r{ui-U2)rduidu2<2tv \r{u)\''+^ du. 
Jo Jo 

Hence, we get from (2.1) and (4.10) that, for all u > 0, as t — > oo, 

rvt rvt 

t— > oo. 

Now pick / E L'^{y) with Ef{No) = and Hermite rank q{f) = q; by (A. 2), 
for ah t > 0, de'N, 9i,...,9deR and Vd > Vd-i > ■ ■ ■ > vi > vq = 0, 

2 



rvt rvt 

(4.11) Mt{v):=sup / \r{ui-U2)rduidu2 = o{tL{t)), 

m>qJo Jo 



^ E^W fiNu)-^Hg{N^)du 

\k=i -^^fe-i* ^- y 

< E yj^ \r{ui-U2)rduidu2 



l<ii,i2<d k=q-\-l 
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l<ii,i2<:d k=q+l 

where are the quantity defined in (1.4). Note that the sum on the right- 
hand side is finite since / G L'^{i^)- We have already shown that (i) holds for 
the particular choice / = TCq. In view of (4.11) and the chain of inequahties in 
the last display, the proof of assertion (i) is completed by the Cramer- Wold 
device. 

(ii) Let / G L'^{i') with Ef{Nu) = and Hermite rank q{f) = q such that 
there exists R>1 with 

oo I I 

(4.12) E7=^'<°°- 

k=q ^ ^■ 

In view of (i), it remains to show that {2tL{t))~^^'^ Jq f{Nu) is tight in 
C(Mq'), as i — > oo. In view of Billingsley's tightness condition ([6], Theo- 
rem 13.14) and the Cauchy-Schwarz inequality, for xq > it suffices to find 
6 > 0, tQ> and C such that, for all < x < xq and t > to, 

/ rxt \ 2(1+0) 

(4.13) E[{2tL{t)r^/^ f{Nu)duj <Cx'+'. 

Therefore, recall the following inequality as shown by Hariz [15] (cf. proof 
of Theorem l(ii) in [15]): for ah < 6* < 1, t > 0, 

/ /•< \2(l+6») r oo om6»/{l+0)| \/pt s 1/2] 2(1+^ ) 

E{lnN.),u) <[^^g^_^(^H„)r.„) _ . 

Consequently, for all < 9 < 1 and t > 0, we have by stationarity and 
Cauchy-Schwarz inequality, 

2(l+6») / , oo om,6l/(l+e) U | \ 2(1+^*) 



f ft \ Z^l+O) / o^ 



"mi 



m=q 

Thus, it follows from (4.12) that there exist 9 G (0, 1] and a constant C such 
that, for all t > 0, 

E^jj{Nu)duJ <C'{2tL\.\{t)f+'. 
Consequently, for all xq > and < x < xq , we have 

/ , r-^i \2(l+6») /Li ifxntlN 

E{{2tL{t)ry^ f{N^)du) <G'(^^) 

As is slowly varying, it follows from (2.1) that there is indeed a constant 
C G M+ and to > such that (4.13) holds for alH > and < x < xq, giving 
(ii)- 
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4.2. Proof of Corollary 2.1. For H G (0, 1), we set = 2/V{2H + 1) 
[cf. (2.11)]. With the help of ^ih, we define a normal process Ni^ = (A''^^)jgK 
by 

(4.14) 7V^^, = 7^4' f e-^^'-'^ dB^, teR, He{0, 1), 7 > 0, 
and set rH,-y (t) = SA^^t^^^o t e I^- Then we have the following lemma: 

Lemma 4.1. Let 7 > 0, H £ (0, 1) and q£'N. 
(i) If 1/2, then 

rH,,{t) = 7^^-^^*^^-^ + 0{t^^-% t ^ 00, 

whereas Pi/2^{i) = e"'^'*' /or all t G R. 
(ii) 

/ rl, (t)dt = 7^isin2(7rii)- / j^——di. 

Jr TT Jo (1 + ? ) 

(iii) 

(4.15) \ngi^^fi]i^Bi^^^)-ngiN^^,)\ds<oo a.s. 

Proof, (i) For H = 1/2, this is a well-known property of the classical 
Ornstein-Uhlenbeck process driven by the Wiener process. For H ^ 1/2, 
assertion (i) is found in Theorem 2.3 in [12]. 

(ii) The spectral density fn of A'^^, that is, rH^^{t) = /j^ e**^/iy_^(^) d^, has 
the following form ([12], (2.2)): 

(4.16) fH,^{0 = i^''^~'M^Hy-^^—^, (em, He (0,1). 

In particular, /^^(^) d^, is finite if and only if if G (0, 3/4). In this case, the 
identity in (ii) follows from Plancherel's theorem and simple substitutions. 

(iii) Fix g G N and H G (0, 1). Then we find a constant Cg G (0,oo) such 
that iHgix) - nq{y)\ < Cq|x - y|(l + Ixl* + |y|9) for ah x,y G M, where Tig is 
defined in (1.2). Note that almost surely, 

(4.17) l^'mB^^t = N^^t - e-^'N:^,o, t > 0, 
which completes the proof of (iii) as almost surely, 

|W,(7^/iH<.)' - H,(iV^^J| ds 

< Cg\N^^Q\ e-^"(l + 2^+ V^'fj" + '^'\Ko\'') ds < 00. 
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We return to the proof of Corollary 2.1: If i7 G (0,3/4) then (2.7) follows 
from the functional CLT by Hariz [15] and Lemma 4.1. By employing the ma- 
chinery developed by Pipiras and Taqqu ([21], Theorem 3.2), it is straight- 
forward to show that Ni^ satisfies the assumptions imposed by Taqqu in 
[25] for ah H G (1/2, 1), and thus (2.7) is in place for all H G (3/4, 1). RecaU 
that the spectral measure of A'^.^ is absolutely continuous with respect to 
Lebesgue measure [cf. (4.16)], and thus FOUP admits the representation in 
(2.2) (cf. Ibragimov and Linnik [16], Theorem 16.7.2). For H = 3/4, (2.7) is 
thus implied by Theorem 2.1 and Lemma 4.1(i). □ 

4.3. Proof of Theorem 3.1. In order to show Theorem 3.1, we need the 
following modification of Corollary 2.1. 

Lemma 4.2. For H G (0,3/4] let be standard normal random vari- 
able and, for H G (3/4, 1), let be a Rosenblatt distributed random vari- 
able with unit variance. Then, for all H G (0, 1), we have 

(4.18) l'^-(^) ^[(O' - '']^(;ify"), t^oo, 



where Y is a standard normal random variable, independent of X^ , and the 
quantities qh and an and are defined in (2.8) and (2.10) and (2.11), 
respectively. 

Proof. Let H G (0, 1) and recall (2.11). It is clear that (S^)^ ^ 
as i — > oo. Hence, weak convergence of the marginals follows from Corol- 
lary 2.1 and Lemma 4.1(ii) and (iii). It remains to show that the marginals 
are asymptotically independent. This is implied by (4.17) and Lemma 4.1(iv), 
provided we can show that the components of the vector in the following 
display are asymptotically independent, as i ^ oo: 

(4.19) ^ 

where Ni. is defined in (4.14). In view of the stationarity of Ni and Lem- 
ma 4.1 (iii), it suffices to show that, for all 6*1,^2 G ]R, 



/•O 

(4.20) lim E 



t—>oo 



0. 



For n G N, let M„ = £^[e*^^^i'"|^n] where J-n is the u-algebra generated 
by {Nig)s<~n- Then {Mn,J^n)n&H is a (bounded) backward martingale. By 
virtue of the convergence theorem for backward martingales (cf. Hall and 
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Heyde [14], Theorem 2.6), note that lim^^^oo Mn = Eexp{i92Nfo) almost 
surely. 

For 7 > 1 pick G N, < t such that nt^oo and gnit) / gH{nt) = as 
t — > oo . In view of the triangular inequality, conditioning on i this yields 
the following inequality, for t > 1 , 



E 



<2E 



By the choice of n^, gnil) Jn^ T~i2{N^) tends to zero in probability for t ^ oo. 
Hence, the right-hand side in the last display tends to zero as t — > oo, giv- 
ing (4.20). □ 

We return to the proof of Theorem 3.1. Fix H G (0, 1), for 7 > 0, it follows 
from the self-similarity of [ci. (2.5)] that, for the finite dimensional 
distributions we have 



(4.21) «,)oo= (i?,^/,-/^-(-*-^)5fJ,rf.)^^^^7-^(i?ff. 



Next, we represent the vector Tnij) in terms of the random quantities on 
the left-hand side in the last display. For 7 > 1 , set s/^ (7) = jgn (7) and 



H _ ^,2H-1 / 1 / TjH \2 



^7 



K<i)' + ^£Ks)' ds) , < = 7^^ /^'KJ^ ds. 



In view of (2.11) and (4.21), we have 

diag{sHh),SH{l), l)(Vf - fiH\W^ - ii-h\tha{i) - 7)' 

lH \2 



(4.22) 



»(7)/;[«)^-.;,v. (1)4^3^111 



fds 



'0^ 



Note that gnil) and 7""*^ !q{Biu)^ du ^Jj^, as 7 ^ 00, by Lemma 4.2. 
Consequently, we get from Lemma 4.2 that, as 7 — > 00, the left-hand side in 
(4.22) converges jointly in distribution to the random vector 



(4.23) 
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Applying the delta method to 'tp{v,w) = {w~^^'^ ,w^^ jVw"^/"^)' , as 7 — > 00, 

SH (7) [diag (7-^ , 7"''' , 7""^ ' ) {th^ (7) , rH,2 (7) , r//,3 (7) )' " if^'n . /^h )] 

= SH(7)[V'(Vf ,<) - V'lMnSA^i^')] -^HXV(Mi^\A^i^') (J) • 

As we have joint convergence in (4.22) to the vector in (4.23), the proof of 
Theorem 3.1 is completed by elementary calculus. 

4.4. Proof of Theorem 3.2. Fix H G (0, 1). For all 7 < 0, we set 
(4.24) 

rn\ 
Jo 

It follows from (2.5) that the following identity holds for the finite dimen- 
sional distributions: 

iB^,^t)t>o = l7r''(5^i,|,|i)t>o, 7 < 0, G (0, 1). 
Consequently, for 7 < 0, 

diag(|7r(^^+^)/^el^l,|7|-^^-^e2H,|7|(^^-^)/M7|-'eH)r^(7) 
= diag(el^l,e^'^',l,el'^l) 

1 1 Th,i(7) TH,ih)\' 



.Th,2{iY'^ ' ^^,2(7) ' {Th,2{i)Y'^ ' ThAi) . 
Integrating by parts and summing and subtracting terms, for 7 < 0, we have 

(4.25) _ 

+ e'^r_ff,i(7)) 

where we set 



Th,i{i)= e-'B^ e-'B^dtds+l{B^^^y- [B^Y ds 

Jo Js JO 

PCX) roo 

-el^la^l/ e-'B^ds + Bf, e-'B^ds, 7 < 0. 

J\-y\ Jo 

As B^^ /t — > for t — > 00 and H ^ (0, 1), almost surely, we note that Cu = 
supi>o ]^_^Y is finite, almost surely. Consequently, we have the fohowing 
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chain of inequalities [7 < 0] , almost surely, 

/ e'B^ e-^Bl^ dtds <Cjj e'{l + s) e-\l + t)dtds 

Jo Js Jo Js 

fill /"oo 

= Cji (l + s) e-\l + s + t)dtds<C'C]j{l + \-ff), 
Jo Jo 

where C" > is a finite deterministic constant not depending on 7. 

Exploring similar bounds for the other terms in (4.26) yields, almost 
surely, 

\fH,i{l)\<C"me.x{CH,Cjj)il + hf), 

where C" is a finite and deterministic constant not depending on 7. Conse- 
quently, almost surely. 



(4.26) 
Similar, 

(4.27) 



6^7^^,2(7)^0, 



7" 



-00. 



Ip2|7l 







00 \ 2 

e-'B^ds) +Ti^,3(7), 



7<0, 



where e^'^Tn^zil) 0, for 7 —00, almost surely. 
Finally, for all H S (0, 1) and 7 —00, note that 



(4.28) 



V 



e-'B^ds 



I 



Z 



where (y, Z) is standard normal. Indeed, as all corresponding vectors are cen- 
tered Gaussian, it suffices to investigate their covariance matrices. Clearly, 
/o°° ds is Gaussian with a mean zero and variance T{2H + l)/2. Fur- 

thermore, 

j\"'+'{Bf^^ - B^) ds = B^ ds, 7<0. 

The right-hand side tends to e~^Bf ds, almost surely and in mean square 
as 7 — > —00. 

For all 7 G M and H & (0, 1), we have the following identity, 
e(^I^^ e~'B^dt) dJ"' e^+^(i?|^| - ) ds^ 

= e-* e-^i?i?f (5|^| - J ds dt. 
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In view of (2.4), expanding the integrand into powers of s/|7|, t/\^\ and 
(s + i)/|7| up to the second order, this shows that the integrand tends to 
zero for 7^ -oo ands,t eM. Note that \EBl^ {B^^^- B^^^_^)\ < \t\'^^\s\'^" by 
means of the Cauchy-Schwarz inequahty. Thus, the dominated convergence 
theorem is apphcable to the right-hand side in the last display giving (4.28) 
for all H e (0, 1) and 7 -00. In view of (4.25), (4.26) and (4.27), (i) now 
follows from (4.28). 

APPENDIX: THE DIAGRAM FORMULA 

We refer to BoUobas [5] for basic notation and facts of graph theory. The 
sets of vertices and edges of an undirected graph D are denoted by V{D) 
and S{D), respectively. Throughout, we identify S{D) with a subset of the 
set {{v, w} -.VjW G V{D),v / w}. The degree of a vertex v G V{D) is defined 
to be the cardinality of the set {e G £{D) -.vGe}. 

In the sequel, let p, G N. A diagram D with p levels and rows of length q is 
an undirected graph with vertex set V{D) = {1, . . . ,p} x {1, . . . , q}, satisfying 
the following properties: 

(i) U{{h,h),{l2,k2)}e£{D) thenh^h. 

(ii) Each vertex has degree one. 

The (possibly empty) set of all diagrams with p levels and rows of length q 
is denoted by 'D{p, q). For D £ 'D{p, q) and e = {{k, m), {I, n)} £ £{D), set 

(A.l) m{e) =mDie) = inm{k,l}, M{e) = Mnie) = inax{k,l}. 

Now let TCq be the Hermite polynomial of order q in (1.2) and {Xi, . . . ,Xp) 
be a Gaussian vector with EXj = and EXj =1 (1 < j < p). It follows from 
the diagram formula (cf. Breuer and Major [7]) that we have 

(A.2) Ef[n,{X,)= ^ n E[X^ie)XM(e)], 

j=l D&V{p,q) eef (D) 

with the conventions Ylz = 1 a^d = 0- 
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